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The  mechanical  properties  of  light  have  found  widespread  use  in  the  manipulation  of  gas-phase 
atoms  and  ions,  helping  create  new  states  of  matter  and  realize  complex  quantum  interactions. 
The  field  of  cavity-optomechanics  strives  to  scale  this  interaction  to  much  larger,  even  human-sized 
mechanical  objects.  Going  beyond  the  canonical  Fabry-Perot  cavity  with  a  movable  mirror,  here 
we  explore  a  new  paradigm  in  which  multiple  cavity-optomechanical  elements  are  wired  together  to 
form  optomechanical  circuits.  Using  a  pair  of  optomechanical  cavities  coupled  together  via  a  phonon 
waveguide  we  demonstrate  a  tunable  delay  and  filter  for  microwave-over-optical  signal  processing.  In 
addition,  we  realize  a  tight-binding  form  of  mechanical  coupling  between  distant  optomechanical  cav¬ 
ities,  leading  to  direct  phonon  exchange  without  dissipation  in  the  waveguide.  These  measurements 
indicate  the  feasibility  of  phonon-routing  based  information  processing  in  optomechanical  crystal 
circuitry,  and  further,  to  the  possibility  of  realizing  topological  phases  of  photons  and  phonons  in 
optomechanical  cavity  lattices. 


I.  INTRODUCTION 

Microscopic  optical  and  microwave  cavities,  with  their 
wavelength  or  even  sub-wavelength  mode  size,  have  re¬ 
cently  been  used  to  greatly  enhance  the  radiation  pres¬ 
sure  interaction  between  electromagnetic  waves  trapped 
in  the  cavity  and  mechanical  vibrations  of  the  cavity 
walls  pQ.  Technical  advances  in  our  ability  to  fabricate 
structures  of  ever  smaller  size  and  higher  quality  has  led 
to  the  demonstration  of  a  number  of  new  phenomena,  in¬ 
cluding  laser  cooling  of  mechanical  resonators  into  their 
quantum  ground  state  m  ei,  optomechanically-induced 
electromagnetic  transparency  [Ml,  and  microwave  sig¬ 
nal  processing  To  date  these  experiments  have 

largely  involved  single  cavity  elements.  A  new  paradigm, 
enabled  by  chip-scale  cavity-optomechanics,  involves  the 
integration  of  multiple  optomechanical  cavities  together 
to  realize  extended  photon-phonon  excitations  for  classi¬ 
cal  and  quantum  information  processing  applications  0- 
ITTj.  Among  the  different  types  of  optomechanical  sys¬ 
tems  currently  being  studied,  the  device  architecture 
based  on  thin-film  optomechanical  crystals  CMH  is 
a  particularly  promising  integration  platform  given  the 
large  attainable  radiation-pressure  coupling  and  the  abil¬ 
ity  to  create  photonic  and  phononic  band-gap  waveg¬ 
uides. 

In  an  integrated  optomechanical  network  composed  of 
cavities  and  waveguides,  propagating  phonons  can  be 
routed  and  stored  among  different  nodes  as  controlled 
by  driving  laser  fields,  providing  various  information  pro¬ 
cessing  functions  -  such  as  buffering,  delay,  and  filtering  - 
in  a  chip-scale  footprint  due  to  the  slow  velocity  of  acous¬ 
tic  waves  PUCES].  As  an  initial  step  towards  realizing 
more  complex  multi-element  cavity-optomechanical  sys¬ 
tems,  here  we  explore  the  optical  excitation  and  routing 
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of  GHz  phonons  in  a  small  optomechanical  circuit  con¬ 
sisting  of  two  optomechanical  crystal  cavities  connected 
by  a  dispersion-engineered  phonon  waveguide.  Pulsed 
and  continuous- wave  measurements  are  first  used  to  char¬ 
acterize  the  coupling  efficiency  and  phonon  propagation 
properties  of  the  system.  Utilizing  separate  optical  driv¬ 
ing  fields  for  each  optomechanical  cavity,  we  use  this  sys¬ 
tem  to  demonstrate  an  optically-tunable  microwave  de¬ 
lay  line  in  which  microwave-over-optical  signals  sent  in 
one  cavity  port  are  efficiently  converted  into  outgoing 
microwave-over-optical  signals  on  the  other  cavity  port. 


Owing  to  the  continuum  of  modes  supported  by  the 
phonon  waveguide,  a  tight-binding  type  of  interaction 
between  distant  cavities  may  also  be  realized.  Such  ideas 
have  recently  been  explored  utilizing  photonic  waveg¬ 
uides,  both  in  the  optical  domain  with  photonic  crystal 
cavities  |20j  and  in  the  microwave  domain  with  super¬ 
conducting  qubits  [2lj .  Here  we  use  a  phonon  waveg¬ 
uide  to  strongly  couple  two  distant  optomechanical  cav¬ 
ities,  leading  to  an  exchange-type  phonon  routing.  The 
hybridized  mechanical  cavity  modes,  while  spatially  ex¬ 
tended,  nevertheless  have  strong  optomechanical  cou¬ 
pling  with  both  localized  optical  cavity  modes  which  are 
not  hybridized.  Such  a  feature  can  be  used  to  create 
mechanical  transducers  capable  of  bridging  two  sepa¬ 
rated  physical  systems  (22).  such  as  in  quantum  wave¬ 
length  conversion  between  optical  and  superconducting 
microwave  systems  ]2.‘{l-l2fil  and  coupling  between  solid 
state  spins  [271429] .  Scaling  such  strong  coupling  to  arrays 
of  cavity  elements  may  also  enable  novel  optomechani¬ 
cal  metamaterials,  where  in  combination  with  a  spatially 
distributed  phase  pattern  of  the  optical  driving  field  [30] , 
topological  phases  for  photons  and  phonons  are  expected 
to  arise  [33]  [32]  ■ 
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FIG.  1.  GHz  dispersion-engineerable  phonon  waveguide,  a,  A  schematic  diagram  of  the  single  cavity-waveguide  device. 
a  and  b  represent  the  optical  and  mechanical  cavity  modes.  The  mechanical  cavity  mode  couples  to  the  phonon  waveguide  at 
rate  ye  and  has  an  optically  induced  decay  rate  7om-  b,  Scanning  electron  microscope  (SEM)  image  of  a  silicon  optomechanical 
cavity-waveguide  device  based  on  the  scheme  of  a.  c,  Theoretical  band  structure  of  the  phonon  waveguide  employed  in  this 
work.  Green  curve  is  the  waveguide  breathing  mode  and  yellow  curves  are  the  bands  of  the  phononic  shielding.  Inset  on  the 
right  is  the  modal  profile  of  the  waveguide  breathing  mode  corresponding  to  wavevector  k  =  0.8  J.  d,  Optically  transduced 
mechanical  power  spectral  density  (PSD)  of  a  typical  device  showing  a  few  hybridized  cavity-waveguide  resonances,  e,  Group 
velocity  as  a  function  of  frequency  across  the  breathing  mode  waveguide  band.  Green  is  the  theoretical  curve  and  red  dots  are 
the  measured  data. 


II.  PHONON  WAVEGUIDE 
CHARACTERIZATION 

While  integrated  phonon  waveguides  with  frequency 
bands  of  kHz  to  a  few  hundred  MHz  have  been  pre¬ 
viously  studied  J23H2SF  here  we  explore  the  properties 
of  thin-film,  dispersion  engineered  phonon  waveguides 
in  the  GHz  frequency  range.  A  schematic  diagram  and 
SEM  image  of  a  cavity-waveguide  device  used  to  study 
such  phonon  waveguide  properties  are  shown  in  Fig.  Eh 
and  Fig.  [ljj,  respectively.  The  device  is  fabricated  from 
the  220  nrn  silicon  device  layer  of  a  silicon-on-insulator 
wafer  (see  App.  0>  and  consists  of  several  parts:  an  op¬ 
tomechanical  cavity  with  co-localized  optical  and  acous¬ 
tic  resonances,  an  optical  coupling  waveguide,  and  a  long 
phonon  waveguide  section  which  is  end-fire  coupled  to  the 
cavity.  The  optomechanical  cavity  is  realized  in  a  pho¬ 
tonic  crystal  nanobeam  which  supports  an  optical  reso¬ 
nance  at  Ac  ss  1550  nm  in  the  telecom  wavelength  band 
that  is  strongly  coupled  through  the  elasto-optic  radi¬ 
ation  pressure  effect  to  a  “breathing  mode”  mechanical 
resonance  of  the  beam  at  a  frequency  of  uim/ 2tt  «  6  GHz. 
One  end  of  the  nanobeam  is  terminated  by  a  two  di¬ 
mensional  phononic  crystal  mirror  with  a  full  phononic 
bandgap  in  the  relevant  frequency  range  (15) :  the  other 
end  is  connected  to  the  phonon  waveguide  consisting  of 


a  similarly  patterned  central  beam  region  surrounded  by 
phononic  bandgap  shielding  on  either  side.  The  phonon 
waveguide  is  designed  to  have  large  extrinsic  coupling 
7e  to  the  nanobeam  cavity  mode,  leading  to  hybridiza¬ 
tion  of  the  localized  cavity  mode  and  the  nearly-resonant 
phonon  waveguide  modes.  The  optical  coupling  waveg¬ 
uide  is  fabricated  in  the  near-field  of  the  nanobeam  cavity 
(bottom  beam  of  Fig.  E>).  allowing  for  evanescent  cou¬ 
pling  of  laser  light  into  and  out  of  the  cavity.  A  single 
optical  fiber  taper  is  used  to  couple  light  into  the  on- 
chip  coupling  waveguide,  and  a  photonic  crystal  mirror 
is  etched  in  to  the  end  of  the  optical  coupling  waveguide 
so  that  light  coupled  into  the  nanobeam  cavity  can  be 
recollected  by  the  optical  fiber  taper  as  per  Ref.  [36], 

Figure  [TJ-  shows  the  band  structure  of  the  phonon 
waveguide  with  silicon  slab  thickness  d  =  220  nm  and 
lattice  constant  a  =  ao  =  480  nm,  numerically  simu¬ 
lated  using  the  finite-element  method  m-  For  clarity, 
only  the  bands  with  even  z— symmetry  (out-of-plane  di¬ 
rection)  are  shown,  with  the  green  curve  corresponding 
to  the  breathing  mode  waveguide  band.  For  a  lattice 
constant  of  ao  =  480  nm  this  band  has  a  bandwidth  of 
approximately  1  GHz  surrounding  the  6  GHz  frequency 
of  the  localized  acoustic  mode  of  the  nanobeam  cavity. 
Scaling  of  the  planar  waveguide  dimensions  by  scale  fac¬ 
tor  s'  =  a/oo,  where  a  is  the  new  lattice  constant,  re- 
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suits  in  an  approximate  scaling  of  the  waveguide  bands 
by  (sr1-  This  is  a  good  approximation  for  the  lower  ly¬ 
ing  bands  of  thin  structures  in  which  d/a 0  <  1.  As  such, 
here  we  plot  the  bandstructure  in  terms  of  a  (unitless) 
normalized  frequency  s  =  (// 6  GHz)s',  where  /  is  the 
physical  waveguide  band  frequency.  With  the  phonon 
waveguide  terminated  on  one  end  by  the  nanobeam  cav¬ 
ity  and  on  the  other  end  by  cutting  out  a  slot  in  the  center 
beam,  the  continuous  waveguide  band  is  transformed  into 
discrete  waveguide  resonances  with  free  spectral  range 
/fsr  =  vg/( 20>  where  vg  =  \du>/dk\  is  the  group  ve¬ 
locity  of  the  band  and  l  is  the  length  of  the  waveguide. 
For  those  waveguide  modes  spectrally  located  within  a 
bandwidth  ye  of  the  localized  cavity  mode,  optical  exci¬ 
tation  of  the  optomechanical  cavity  can  be  used  to  detect 
phonons  propagating  in  the  waveguide  section  of  the  de¬ 
vice. 

Optical  excitation  and  detection  of  the  waveguide 
phonon  resonances  are  performed  using  a  tunable  nar¬ 
rowband  laser  with  frequency  ( ujp )  blue-detuned  from  the 
optical  cavity  frequency  (wc)  by  the  mechanical  cavity 
frequency  (wm),  A  =  ujp  —  ujc  «  wm.  Efficient  Stokes 
scattering  of  the  laser  pump  field  occurs  for  the  mechan¬ 
ical  breathing  mode  resonances  of  the  cavity-waveguide 
system  which  lie  within  an  optical  cavity  linewidth  (re)  of 
the  bare  mechanical  frequency  uim.  For  the  blue-detuned 
sideband  pumping  employed  throughout  this  work,  the 
mechanical  modes  also  experience  amplification  by  the 
optical  pump  wave  corresponding  to  a  negative  (phonon) 
decay  rate  of  7om  =  — Agaric/ n,  where  go  is  the  vacuum 
optomechanical  coupling  rate  characterizing  the  strength 
of  the  radiation  pressure  interaction  [5] ,  nc  is  intra-cavity 
photon  number  determined  by  the  optical  pump  power, 
and  re  is  optical  cavity  damping  rate. 

All  experiments  are  performed  at  room  temperature 
and  atmospheric  pressure.  The  reflected  optical  signal 
from  the  cavity  is  detected  on  a  photodetector,  with 
the  beating  of  the  pump  wave  and  resulting  Stokes 
waves  yielding  a  photocurrent  spectrum  containing  the 
(mostly)  thermal  motion  of  the  mechanical  resonances  of 
the  cavity- waveguide  system.  Figure  E*  shows  a  typ¬ 
ical  microwave  photocurrent  spectrum  of  a  fabricated 
device  with  a  waveguide  length  l  =  145  /jm,  cavity-to- 
waveguide  coupling  rate  7e/27r  =  2.4  MHz,  and  loaded 
optical  cavity  linewidth  re/27 r  =  0.8  GHz.  A  series  of  hy¬ 
bridized  modes  with  free-spectral  range  /fsr  =  4  MHz 
are  observed  around  the  main  cavity  resonance  peak  at 
wm/27r  =  6.017  GHz,  from  which  one  can  infer  a  waveg¬ 
uide  group  velocity  of  vg  =  1160  m/s.  In  order  to  trace 
out  the  group  velocity  evolution  along  the  entire  breath¬ 
ing  mode  waveguide  band  a  set  of  devices  is  fabricated 
with  the  planar  dimensions  of  the  waveguide  scaled  by 
different  scaling  factors  s'  =  a/ao ■  The  nanobeam  cav¬ 
ity  parameters  are  held  fixed.  In  this  way  the  waveguide 
band  can  be  shifted  in  frequency  relative  to  the  nomi¬ 
nally  fixed  cavity,  allowing  us  to  effectively  scan  across 
the  whole  waveguide  band  using  the  cavity  as  a  frequency 
(and  spatial  mode)  filter.  Figure  [jj  shows  the  measured 


group  velocity  evolution  along  the  waveguide  band  as  in¬ 
ferred  from  the  measured  /fsr  of  the  waveguide  reso¬ 
nances  versus  normalized  frequency  s.  The  green  curve 
corresponds  to  the  theoretically  calculated  group  velocity 
curve.  We  note  that  in  the  lower  frequency  portion  of  the 
band  there  exists  two  fc-points  for  each  frequency,  result¬ 
ing  in  two  group  velocity  curves,  a  feature  also  seen  in  the 
measured  data.  In  analogy  to  the  slow  light  in  photonic 
crystals  m,  here  we  measure  slow  phonon  group  veloc¬ 
ities  [35]  near  both  band  edges  as  well  as  at  the  band 
minimum  occurring  near  k  «  0.5(7r/a).  The  minimum 
measured  phonon  group  velocity  is  as  small  as  150  m/s, 
a  reduction  by  a  factor  of  40  compared  to  the  group  ve¬ 
locity  of  transverse  acoustic  waves  in  bulk  silicon  at  room 
temperature. 


III.  OPTOMECHANICAL  MICROWAVE  DELAY 
LINE  AND  FILTER 

We  next  study  the  utility  of  cavity-waveguide  optome¬ 
chanical  systems  for  coherent  microwave-over-optical  sig¬ 
nal  processing.  Similar  signal  processing  based  on  stimu¬ 
lated  Brillouin  scattering  has  recently  been  investigated 
in  photonic  nanoscale  waveguides  [88Tt40] .  Here  we  em¬ 
ploy  cavities  as  localized  sites  for  enhancing  the  interac¬ 
tion  of  photons  and  phonons,  and  phonon  waveguides  for 
the  routing  and  filtering  of  microwave  signals.  The  ba¬ 
sic  circuit,  shown  schematically  in  Fig.[2ji,  consists  of  two 
optomechanical  cavities,  labeled  L  and  R ,  which  are  con¬ 
nected  by  a  common  phonon  waveguide.  The  microwave 
signal  encoded  in  the  single  sideband  of  an  optical  car¬ 
rier  with  wavelength  Al  ,  is  fed  into  cavity  L  which  coher¬ 
ently  drives  the  localized  cavity  breathing  mode  through 
the  radiation-pressure  force.  The  excited  phonons  then 
couple  into  the  phonon  waveguide  and  propagate  to  cav¬ 
ity  R,  providing  the  desired  delay  and  filter  functions  for 
the  microwave  signal.  In  cavity  R ,  the  arriving  phonons 
modulate  the  internal  optical  field,  re-encoding  the  mi¬ 
crowave  signal  on  an  outgoing  carrier  at  wavelength  A«. 

A  schematic  of  the  experimental  set-up  used  to  mea¬ 
sure  the  two-cavity  system  is  shown  in  Fig.  [2jo.  SEM 
images  of  the  cavity  and  waveguide  elements  as  they 
are  arranged  in  the  optomechanical  circuit  are  shown  in 
Fig.  An  adiabatic  fiber-to-chip  coupling  waveguide 
followed  by  a  ’Y’-branch  splitter  is  used  to  couple  laser 
light  into  and  out  of  each  of  the  cavity  devices  using  a 
single  optical  fiber  taper.  The  two  nanobeams  are  de¬ 
signed  to  have  similar  mechanical  cavity  frequencies  but 
different  optical  resonant  frequencies  so  that  they  may 
be  simultaneously  and  independently  probed  using  two 
separate  laser  sources  (see  App.  0' 

A  typical  mechanical  spectrum  of  the  double-cavity 
system  with  phonon  waveguide  length  l  =  43  p,m  is 
shown  in  Fig.  11-  For  this  device  the  two  funda¬ 
mental  optical  cavity  modes  Ol(,r)  have  wavelength 
^L(R)  =  1539(1559)  nm  and  linewidth  Kl{r)  =  27t  x 
1.15(0.80)  GHz.  The  blue  curve  corresponds  to  blue- 
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FIG.  2.  Microwave  signal  processing  using  optomechanical  cavity-waveguide  system,  a,  An  optomechanical  scheme 
for  microwave  signal  processing  relying  on  phonon  routing,  b,  Schematic  of  the  experimental  set-up.  VOA,  variable  optical 
attenuator;  FPC,  fiber  polarizer  controller;  PG,  pulse  generator;  EOM,  elector-optic  modulator;  A-(D)mux,  wavelength  (de¬ 
multiplexer;  PD,  photodetector;  EDFA,  erbium  doped  fiber  amplifier;  SA,  spectrum  analyzer;  NA,  network  analyzer;  BPF, 
bandpass  filter;  OS,  oscilloscope,  c,  SEM  image  of  the  double-cavity  device,  showing  the  ’Y’-branch  optical  splitter  which  feeds 
both  left  and  right  optical  cavities  and  the  intermediary  mechanical  waveguide  that  couples  the  left  and  right  cavities  together. 
The  mechanical  spectrum  of  the  device  as  measured  from  cavity  L  (blue  curve)  and  R  (orange  curve)  is  shown  in  the  top  inset. 


detuned  sideband  pumping  of  cavity  L  (Ar/2tt  =  6  GHz; 
Ppl  =  180  /zW;  nCL  =  1200)  and  the  laser  resonant 
with  cavity  R  turned  off.  The  orange  curve  is  for  similar 
pumping  of  cavity  R  (Ar/2tt  =  6  GHz;  Ppr  =  50  /z W; 
nCR  =  360).  The  two  spectra  are  seen  to  share  a  number 
of  common  resonances,  a  result  of  the  connecting  phonon 
waveguide.  We  label  the  mechanical  resonances  as  Lj 
( Rk ),  with  the  index  j  ( k )  indicating  the  relative  optical 
coupling  strength  of  the  resonance  to  cavity  L  ( R ).  The 
most  intense  mechanical  resonance  peaks  in  the  cavity 
spectra,  L\  (=  R4)  and  Ri,  are  found  to  have  vacuum 
optomechanical  coupling  rates  of  goLi/27r  =  0.85  MHz 
to  optical  mode  Ol  and  g0Rx /2tt  =  1.39  MHz  to  op¬ 
tical  mode  Or,  respectively,  determined  from  the  vari¬ 
ation  of  the  mechanical  linewidth  with  pump  power. 
The  intrinsic  Q-factor  of  the  measured  mechanical  res¬ 
onances  are  all  approximately  1500  as  determined  from 
the  resonance  peak  linewidths  at  low  optical  probe  power 
(7i/27r  «  3.7  MHz). 

We  next  study  the  propagation  of  short  phonon  pulses 
within  the  double-cavity  system.  In  this  experiment 
the  two  cavities  are  simultaneously  being  probed  with 
two  separate  continuous  wave  lasers  (Ppl  =  80  /zW, 


nCL  =  530;  Ppr  =  50  /zW,  nCR  =  360),  each  tuned  to 
the  blue  sideband  of  their  respective  cavity  resonance. 
A  20  ns  microwave  pulse  is  used  to  phase  modulate  the 
laser  light  input  to  cavity  L  at  a  frequency  of  5.952  GHz 
resonant  with  the  L\  mechanical  resonance,  exciting  a 
coherent  phonon  population  in  cavity  L  (see  App.  0- 
Figure  [3|i  shows  the  optically-transduced  signal  in  both 
cavities  after  the  initial  excitation  pulse.  For  pulses  with 
bandwidth  larger  than  the  free-spectral  range  (FSR)  of 
the  waveguide  modes  connecting  the  two  cavities,  the 
delay  associated  with  phonons  traversing  the  waveguide 
can  be  resolved  and  is  given  by,  r  =  1/2/fsr-  Multiple 
bounces  of  the  initial  phonon  pulse  are  seen  in  the  sig¬ 
nal  from  each  cavity,  with  a  period  between  bounces  of 
58  ns  (corresponding  to  r  =  29  ns)  consistent  with  the 
measured  FSR  of  17  MHz.  The  pulse  propagation  can 
be  simulated  with  coupled  mode  equations  (see  App. 
and  are  shown  as  dashed  lines  in  Fig.  UK  From  a  fit 
of  the  model  to  the  measured  pulsed  response,  the  ex¬ 
trinsic  mechanical  cavity-to-waveguide  coupling  rate  is 
estimated  to  be  7eL(R)/27r  =  35(26)  MHz  for  the  L  ( R ) 
cavity,  corresponding  to  an  internal  efficiency  of  phonon 
transmission  from  cavity  L  to  R  of  89%  (67%)  excluding 
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FIG.  3.  Pulsed  and  CW  signal  propagation,  a,  A  short 
phonon  pulse  excited  by  a  20  ns  microwave  pulse  of  frequency 
5.952  GHz  in  cavity  L  is  detected  in  the  two  optomechan¬ 
ical  cavities  as  it  bounces  between  them.  Dashed  lines  are 
simulated  results,  b,  CW  microwave  propagation  from  cavity 
R  to  cavity  L  is  characterized  by  S'-matrix  component  Srl 
as  measured  by  a  network  analyzer  (see  App.[A|.  Blue  curve 
is  \Srl\2  (in  dB)  and  green  curve  is  the  phase  of  Srl ■  c, 
Microwave  delay  at  peak  of  f?i  resonance  (5.991  GHz)  ver¬ 
sus  pump  power  in  cavity  R  ( Ppr )  with  fixed  pump  power 
in  cavity  L  (Ppl  =  0.2  mW)  and  microwave  signal  power 
P3  =  —20  dBm.  d,  |Sr.l|2  spectrum  around  Ri  resonance  cor¬ 
responding  to  peak  delay  of  13.3  [is.  Red  curve  is  a  Lorentzian 
fit. 


(including)  phonon  waveguide  decay. 

Continuous- wave  (CW)  measurements  of  the  same  de¬ 
vice  are  presented  in  Fig.  §>  in  which  the  frequency  (u>) 
of  the  phase  modulation  applied  to  the  laser  light  sent 
into  cavity  L  is  varied.  The  displayed  S'— matrix  com¬ 
ponent  Srl  (see  App.  E  is  for  an  average  optical  pump 
power  of  PpL(R)  =  80  /iW  for  both  cavities.  Signal  de¬ 
lay  in  this  case  can  be  inferred  from  the  phase  <f>  of  Srl , 
r  =  —  where  r  now  includes  both  the  delay  asso¬ 
ciated  with  propagation  through  the  phonon  waveguide 
and  the  resonant  storage  time  in  each  cavity.  Unlike  pre¬ 
viously  studied  optomechanically-induced  transparency 
in  a  single  cavity  SHI!,  where  delay  is  associated  with 


reduced  transmission  because  of  the  destructive  interfer¬ 
ence  between  direct  optical  transmission  and  phonon  me¬ 
diated  optical  transmission,  here  the  delay  purely  relies 
on  phonons  as  no  light  directly  propagates  through  the 
waveguide  and  one  can  simultaneously  realize  high  trans¬ 
mission  efficiency  and  large  delay.  For  blue-detuned  side¬ 
band  pumping,  which  results  in  parametric  amplification 
of  the  phonon  signal,  the  peak  optical  transmission  gain 
between  cavities  for  the  phonon-waveguide  mediated  pro¬ 
cess  can  be  greater  than  unity,  and  is  given  for  excitation 
of  mechanical  resonance  Lj  (=  R *.)  by, 


G 


j,k,  max  — 


4  CLjCRk 


(i  -cLj-cRkr 


(1) 


where  CLj(Rk)  =  |70M,i3(Hfc) |/7i  is  the  cooperativity  be¬ 
tween  the  Lj  (=  Rk)  mechanical  resonance  and  the  Ol(R) 
optical  cavity  mode.  This  is  the  internal  gain,  and  ne¬ 
glects  optical  loss  outside  of  the  cavities  and  photodetec¬ 
tor  inefficiency. 

Figure  &  shows  the  measured  delay  for  a  5.991  GHz 
microwave  signal  resonant  with  the  R\  mode  as  a  func¬ 
tion  of  the  optical  pump  power  sent  into  cavity  R  for 
fixed  pump  power  Ppl  =  0.2  mW  in  cavity  L.  The  delay 
is  tunable  with  optical  pump  power,  reaching  a  value  as 
large  as  r  =  13.3  / is  for  PpR  =  0.21  mW.  Figure  [3Ji  dis¬ 
plays  the  measured  S— matrix  amplitude  spectrum  cor¬ 
responding  to  the  longest  measured  delay,  showing  a  nar¬ 
row  3-dB  bandwidth  of  17  kHz  and  a  peak  internal  gain 
of  3  dB.  Viewed  as  a  tunable  bandwidth  microwave  fil¬ 
ter  mmm,  this  device  realizes  an  out  of  band  rejec¬ 
tion  of  up  to  70  dB  and  a  thermal-noise-limited  sensitiv¬ 
ity  of  —30  dBm  referred  to  the  microwave  input  of  the 
electro-optic  phase  modulator  (see  App.  0- 


IV.  WAVEGUIDE-MEDIATED  DISTANT 
MECHANICAL  COUPLING 


In  addition  to  guiding  the  propagation  of  acoustic 
waves,  phonon  waveguides  may  also  be  used  to  medi¬ 
ate  coupling  between  distant  optomechanical  cavities  via 
virtual  phonons,  i.e.  without  energy  distribution  in  the 
mediating  waveguide  modes.  Such  a  scenario  can  be 
used  to  realize  highly  efficient  signal  conversion  between 
cavity  elements,  or  in  the  case  of  mechanical  sensing, 
may  be  employed  to  realize  new  sensing  modalities  in 
which  read-out  and  sensor  are  separated  spatially.  As 
recently  demonstrated  for  optical  photons  in  photonic 
crystals  EDI,  when  the  cavity-to- waveguide  coupling  rate 
is  substantially  smaller  than  the  FSR  of  the  waveguide 
modes,  7e/27r  <C  /fsr>  the  excitation  of  real  waveg¬ 
uide  phonons  is  suppressed  and  only  second-order  virtual 
transitions  between  cavity  and  waveguide  modes  con¬ 
tribute  to  the  inter-cavity  coupling.  For  nearly  degen¬ 
erate  bare  cavity  modes  the  waveguide  mediated  inter¬ 
cavity  coupling  strength  scales  as  ye  =  y/'yeL”/eR,  and  for 
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FIG.  4.  Waveguide-mediated  distant  mechanical  cou¬ 
pling.  a  Mechanical  spectrum  of  a  device  with  two  distantly- 
coupled  cavities  as  measured  from  cavity  L  (blue)  and  R  (or¬ 
ange).  Modes  C+  and  C_  are  the  (approximately)  even  and 
odd  hybridized  mechanical  cavity  modes  with  the  simulated 
modal  profile  shown  in  the  inset.  The  weakly  transduced  W\ 
and  W2  resonances  are  waveguide  modes,  b  Measured  phonon 
Rabi  oscillation  between  the  two  distant  cavities. 


7e  >  7i+70M  such  a  configuration  can  provide  both  high 
internal  (i.e.  phonon  routing  between  cavities)  and  exter¬ 
nal  (i.e.  phonon-photon  exchange  in  the  cavity)  efficiency 
for  the  microwave  filter  demonstrated  in  the  previous  sec¬ 
tion. 

Figure  [4^,  shows  the  measured  (thermal)  mechanical 
spectra  of  a  device  with  two  optomechanical  cavities  cou¬ 
pled  together  by  a  30  fim  long  phonon  waveguide  which 
is  designed  to  have  7e/2 /pi  -C  /fsr-  This  device  is  simi¬ 
lar  to  that  in  Fig.  with  one  key  difference:  the  central 
holes  of  the  waveguide  have  been  removed,  reducing  the 
cavity-to-waveguide  coupling.  An  optically-transduced 
spectrum  measured  through  both  the  left  (L)  and  right 
( R )  optical  cavities  is  shown.  The  two  hybridized  me¬ 
chanical  cavity  modes,  labeled  C+  and  CL,  have  a  fre¬ 
quency  splitting  A 4 _ / 27 r  =  4  MHz  and  are  separated 

from  the  nearest  waveguide  mode  by  <5  =  17  MHz.  Both 
of  the  C+  and  CL  mechanical  resonances  are  measured 
to  have  large  optomechanical  couplings  of  over  700  kHz 
with  both  optical  cavity  modes,  indicating  strong  mixing 
of  the  L  and  R  mechanical  cavity  modes.  The  weakly 
transduced  waveguide  modes  are  measured  to  have  a  FSR 
of  /fsr  =  54  MHz. 


From  the  matrix  of  optomechanical  rates  between  the 
hybridized  mechanical  modes  and  the  (uncoupled)  op¬ 
tical  cavity  modes  we  determine  a  waveguide-mediated 
coupling  of  V/2ir  =  1.94  ±  0.06  MHz  and  a  bare  mechan¬ 
ical  mode  splitting  of  Alr/2it  =  0.98  ±  0.48  MHz  (see 
App.  |E|.  The  coupling  rate  between  nearly  degenerate 
L  and  R  mechanical  cavity  modes  can  be  shown  to  be 
given  to  lowest  order  by  V  =  7e/ sin(7r<5//FSR),  yield¬ 
ing  an  average  cavity-to-waveguide  coupling  of  7e/27r  ss 
1.62  MHz,  consistent  with  the  numerically-simulated  de¬ 
sign  value  for  this  device.  Despite  the  distant  separation 
between  the  two  cavities,  phonons  can  nevertheless  tun¬ 
nel  between  them  via  virtual  excitation  of  the  waveguide 
modes.  To  show  this,  a  50  ns  pulse  of  a  5.95  GHz  mi¬ 
crowave  modulation  signal  is  applied  to  the  optical  laser 
beam  input  to  optical  cavity  L,  locally  exciting  a  coherent 
phonon  population  in  the  left  cavity.  Following  the  pulse, 
Fig.  [4)3  shows  the  resulting  oscillation  of  the  phonon  pop¬ 
ulation  between  the  two  cavities.  The  measured  period 
of  oscillation  is  250  ns,  corresponding  to  a  Rabi  frequency 
of  4  MHz  commensurate  with  the  hybridized  mode  split¬ 
ting  A_| _ /2ir.  This  result  should  be  contrasted  with  the 

previously  demonstrated  phonon  pulse  bouncing  shown 
in  Fig.  [3j  where  the  period  is  set  by  the  free  spectral 
range  of  the  waveguide  modes. 


V.  DISCUSSION 

The  devices  presented  in  this  work,  involving  pairs  of 
optomechanical  cavity  elements  connected  together  by  a 
phononic  waveguide,  demonstrate  the  feasibility  of  creat¬ 
ing  small  optomechanical  crystal  circuits  for  microwave- 
over-optical  signal  processing  and  novel  sensing  modal¬ 
ities.  Advantages  of  the  optomechanical  approach  are 
myriad.  For  instance,  the  use  of  radiation  pressure  as  op¬ 
posed  to  the  piezo-electric  effect  to  couple  to  microwave 
phonons  allows  one  to  work  with  a  wider  variety  of  ma¬ 
terials,  such  as  silicon  in  the  devices  studied  here,  and  to 
realize  high  efficiency  of  signal  conversion.  Purely  elec¬ 
tromechanical  nanoscale  devices  j43;^  04] ,  by  comparison, 
suffer  from  capacitance  and  impedance  mismatch  which 
limits  efficiency.  Implementation  of  two-dimensional  op¬ 
tomechanical  crystals,  as  has  been  recently  demonstrated 
in  Ref.  [45],  can  also  be  used  to  create  multi-tap  phonon 
waveguide  structures  capable  of  realizing  diverse  spectral 
filter  shapes  and  dispersion  control  without  the  need  for 
additional  electrical  or  photonic  hardware. 

Further  scaling  of  optomechanical  circuits,  to  arrays 
of  a  large  number  of  cavity  elements,  is  also  interest¬ 
ing.  A  straightforward  one-dimensional  extension  of  the 
currently  studied  structures  to  a  chain  of  mechanically- 
coupled  optomechanical  crystal  cavities,  for  instance,  can 
be  used  to  study  mechanical  collective  dynamics  IS]. 
An  even  richer  set  of  physical  phenomena  and  device 
functionality,  however,  may  be  realized  in  optomechan¬ 
ical  circuits  in  which  both  phononic  and  photonic  con¬ 
nectivity  of  cavity  elements  are  employed.  For  exam- 
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pie,  building  on  the  tight-binding  mechanical  coupling 
of  cavities  demonstrated  in  this  work,  an  acoustic  and 
optical  isolator  may  be  realized  using  a  phase-correlated 
optical  pumping  scheme  in  which  the  cavities  are  also 
coupled  together  by  a  photonic  waveguide  [TUI  @7].  In 
this  scheme  [48],  analyzed  in  App.  [Fj  the  interference 
between  optical  and  mechanical  hopping  amplitudes  be¬ 
tween  the  cavities  picks  up  the  phase  difference  between 
the  laser  pump  beams  at  each  cavity  site,  giving  rise  to 
the  required  non-reciprocal  behavior.  With  the  help  of 
post-processing  methods,  such  as  scanning  probe  nano¬ 
oxidation  lithography  [43]  which  has  been  successfully 
employed  to  tune  photonic  crystal  cavities  ED],  the  task 
of  aligning  both  mechanical  and  optical  resonances  of  the 
optomechanical  cavities  could  be  fulfilled.  This  would  al¬ 
low,  at  least  in  moderate  lattice  sizes,  the  study  of  topo¬ 


logical  phases  and  effective  magnetic  fields  for  photons 
and  phonons  [30U32], 
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Appendix  A:  Methods 


3.  Pulse  and  CW  measurements 


1.  Fabrication 

The  devices  are  fabricated  using  a  silicon-on-insulator 
wafer  with  silicon  device  layer  thickness  of  220  nm  and 
buried-oxide  layer  thickness  of  2  /mi.  The  device  geom¬ 
etry  is  defined  by  electron  beam  lithography  followed  by 
inductively-coupled-plasma  reactive  ion  etching  to  trans¬ 
fer  the  pattern  through  the  220-nm  silicon  device  layer. 
The  devices  are  then  undercut  using  an  HFrH^O  solu¬ 
tion  to  remove  the  buried  oxide  layer,  and  cleaned  using 
a  piranha  etch. 


2.  Experimental  Set-Up 

Two  lasers  are  used  to  probe  the  two  nanobeam 
cavities  individually,  one  of  which  can  be  encoded  with 
microwave  signals  by  an  electro-optic  phase  modula¬ 
tor.  The  output  of  the  two  lasers  are  combined  in  a 
wavelength  multiplexer  and  coupled  into  the  device 
through  a  tapered  optical  fiber.  The  on-chip  optical 
waveguide  is  split  into  two  paths  to  transmit  light 
into  the  two  nanobeam  cavities  through  side-coupled 
optical  couplers.  The  reflected  light  from  the  cavities 
is  re-collected  by  the  optical  fiber  and  passed  through 
an  optical  circulator  to  separate  it  from  the  input  light. 
A  fiber-based  wavelength  de-multiplexer  with  20  nm 
channel  spacing  is  then  used  to  separate  the  two  cavity 
signals.  Finally,  the  reflected  light  is  amplified  by  an 
erbium-doped  fiber  amplifier  and  sent  to  a  high  speed 
photodetector  (12  GHz  bandwidth).  The  electrical 
signal  from  the  photodetector  is  measured  either  with 
a  spectrum  analyzer,  network  analyzer,  or  high  speed 
oscilloscope. 


For  the  phonon  pulse  measurement,  a  pulse  genera¬ 
tor  is  used  to  gate  the  microwave  source  which  drives 
the  electro-optic  modulator  with  a  repetition  rate  of 
20  kHz.  The  photodetected  optical  signal  is  mixed  down 
to  300  MHz  with  a  second  microwave  source  which  is 
locked  to  the  one  driving  the  electro-optic  modulator. 
The  down-converted  signal  is  filtered,  amplified  and  av¬ 
eraged  104  times  on  an  oscilloscope  with  sampling  rate  of 
10  GSa/s.  For  the  CW  measurements  a  network  analyzer 
is  employed,  with  the  microwave  output  of  the  analyzer 
(port  1)  connected  to  the  electro-optic  modulator  used 
to  phase  modulate  the  light  sent  into  cavity  L,  and  the 
photodetected  optical  signal  reflected  from  the  cavity  R 
is  connected  to  port  2  of  the  analyzer.  We  define  the  S- 
matrix  component  Srl  =  <§21,  where  §21  has  been  nor¬ 
malized  to  remove  the  effects  of  optical  insertion  loss  into 
the  devices,  optical  pre-amplification,  and  photodetector 
gain. 

Appendix  B:  Simulation  of  phonon  pulse 
propagation 

In  this  section,  we  show  propagation  and  bouncing  of 
phonon  pulses  in  the  cavity- waveguide  system  (Fig.  3a) 
can  be  well  simulated  by  a  group  of  coupled  mode  equa¬ 
tions  using  input-output  formalism.  The  dynamics  cap¬ 
tured  by  the  coupled  mode  equations  is  a  phonon  pulse 
travelling  in  a  waveguide  terminated  by  two  cavities  with 
bare  mechanical  frequency  uJmL,R  and  waveguide  cou¬ 
pling  rate  7 eL,R-  We  approximate  u)mL,R  to  be  the  fre¬ 
quency  of  cavity-dominated  modes  L\  and  R\  in  the  sim¬ 
ulation.  Since  the  response  time  of  the  optical  cavity  is 
much  shorter  than  that  of  the  mechanical  cavity,  we  can 
exclude  the  dynamics  of  optical  modes  from  these  equa¬ 
tions.  Thus,  the  coupled  mode  equations  can  be  written 
as  follows, 


^  =  -(iumL  +  -  «/oLaOL<4Le-^0(f)©(T  -  t)  +  ^7 7Ebin,L(t),  (Bl) 

-  -(iUmR  +  t  +  V^binAt),  (B2) 

bin,L(t)  =  e-al(AYAbR(t  -  tw)  -  bmM{t  ~  tw)),  (B3) 

fcin.fi(t)  =  e~al(Af 7EbL(t  ~  tw)  -  bintL(t  ~  tw )),  (B4) 


where  «ol  and  o+l  are  the  amplitudes  of  optical  pump 
and  its  red  sideband  in  the  left  cavity,  r  is  the  duration  of 
excitation  pulse,  ujs  is  the  frequency  of  pulse,  0(t)  is  the 
Heaviside  step  function,  7  is  the  effective  decay  rate  of 
the  excited  mechanical  mode,  a  «  j/vg  is  the  waveguide 
loss  rate,  and  tw  =  1/(2/FSr)  -  l/(7ei  +7)  -  V (leR  +7) 
is  the  single  trip  time  the  pulse  spent  in  the  waveguide. 


From  the  mechanical  spectrum  we  find  7  =  27r  x  2.1 
MHz  for  L\  mode  (the  main  coherently- driven  mode) 
during  the  pulse  measurement;  and  by  fitting  the  pulse 
tails  detected  in  each  cavity  we  find  7 eL  =  2nx  34.7  MHz 
and  7 cr  =  27 r  x  25.5  MHz.  Using  these  parameters, 
and  \bR\  can  be  numerically  calculated  from  the  coupled 
mode  equations  and  the  proportional  voltage  signals  are 
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TABLE  I.  Mechanical  mode  parameters 


g0/2n  (MHz) 

7i/27T  (MHz) 

7e/27T  (MHz) 

L\ 

0.85 

3.7 

34.7 

0.75 

3.7 

^3 

0.63 

3.9 

Ri 

1.39 

3.6 

25.5 

shown  in  Fig.  3a.  The  simulated  result  captures  the  main 
features  of  the  measured  pulse  data.  In  particular,  the 
pulse  splitting  observed  from  cavity  R  is  due  to  the  fact 
that  the  pulse  frequency  is  not  in  resonant  with  cavity  R 
and  thus  experiences  destructive  interference  inside  this 
cavity. 

The  phonon  transfer  efficiency  from  cavity  L  to  cavity 
R  is  about  e~'Y^2^FSB-'>  ~  67%.  The  phonon  transfer  ef¬ 
ficiency  from  cavity  to  waveguide  for  cavity  L  and  R  is 


7eL(i?)/(7eL(K)  +  7)  ~  94%(92%). 

We  summarize  the  measured  mechanical  mode  param¬ 
eters  of  the  device  in  Fig.  3  in  Table  [TJ  where  g0  for  the 
Lj(Rk)  modes  is  with  respective  to  optical  cavity 

inodes. 


Appendix  C:  Derivation  of  the  microwave  S'  matrix 
for  the  optomechanical  cavity-waveguide  system 

Here,  we  derive  the  S— matrix  for  a  microwave  signal 
traversing  the  optomechanical  cavity-waveguide  system. 
We  assume  the  mechanical  amplitude  is  small  such  that 
only  the  first-order  optical  sideband  needs  to  be  consid¬ 
ered.  In  the  next  section,  we  will  verify  the  small  me¬ 
chanical  amplitude  assumption. 

The  Hamiltonian  of  the  system  under  continuous  wave 
operation  involves  two  optical  cavity  modes  cll,r  with 
frequency  u}crr  parametrically  coupled  to  a  common  me¬ 
chanical  mode  b , 


H  =  ^2  hucka\ak  +  huimtfb  +  ^  %ofc(^  +  b)a\ak  +  ^  ih^K^kapke  lUpkt(ak  -  a\), 


(Cl) 


k—L,R 


k=L,R 


k=L,R 


where  the  last  term  is  the  pumps  of  the  two  optical  cav¬ 
ities.  For  simplicity  we  assume  the  pumping  lasers  are 
blue-detuned  from  the  cavity  resonances  which  is  true 
for  all  of  our  experiments.  Suppose  the  pumping  laser 
for  cavity  L  is  modulated  at  frequency  w  by  a  microwave 
signal,  then  the  operators  of  the  system  can  be  decom¬ 
posed  into  carriers  and  sidebands, 


ak  =  a0ke-iu^  +  a+ke"4*"-*-")* ,  b  =  /Ue"^,  (C2) 

where  we  only  keep  the  red  sideband  of  the  pumping 
lasers  because  of  rotating  wave  approximation,  given  the 
sideband  resolved  condition  wm  nk  of  our  device.  Sup¬ 
pose  the  pumps  are  strong  enough  such  that  the  car¬ 
rier  operators  can  be  treated  as  static  variables,  then  the 
equations  of  motion  of  the  system  can  be  derived  after 


substituting  Eq.  C2  into  Eq.  Cl 


iwa+k  —  ( iAk  „  )o_|_fc  'i'90k&0k{3—  -\/ RekCX m,k> 

(C3) 

-iu/3-  =  -{ium  +  y)/3-  -  igokUokaXk,  (C4) 


where  Ak  =  u}pk  —uick  ~  cum.  Solving  Eq.  C3  and  Eq.  C4 
in  the  frequency  range  | to  —  um\  -C  Kk ,  we  obtain 


P-  = 


WOLy/HeL-^aoL 


i(Um  ~  V)  +  f  ~  Ej 


29nJ“0fcl2 


ah 


(C5) 


^out,i?  —  \/  ^eR^+R 

^SOLgoRy/KeLKeR./ {K'LKr)(XqL(XoR 


i(ujm  -  OJ)  +  f  -  J2k 


29nfe  la0fc  I2 


From  Eq.  C6  peak  optical  gain  at  o j  =  ojm  is 

12  4  CLCR 


C,  v 


|^in,L  |2  (1  —  Cl  —  Cr)2' 


(C6) 

^in,L  • 

(C7) 


where  Cr^r)  =  |7oM,L(fl) \/li  is  the  coorperativity  of  me¬ 
chanical  mode  b  with  optical  modes 


Using  the  result  of  Eq.  C6  the  microwave  signal  trans¬ 
fer  S'— matrix  can  be  derived 


Srl 

Una.i 


(C8) 


Una, out 

7oA7oi?.Ge GeDFA  (J'b^OcR^m/ ^ .  R 

(2Ur/7T)( 

C^in,  l/(*W 

m 

VoLVoRGeGEDFA  ^goLgoRbuJcR^m/ (KLKR)\aoL\2\oiOR\2 


2UW/7T 


i{um  -  ui)  +  ^ 


29ofc  l“0fc  I2 


where  Una, out  and  Una, in  are  the  output  and  detected 
electrical  voltage  of  the  network  analyzer  respectively, 
VoL,r  is  the  optical  loss  of  the  input  and  output  ports  of 
the  device  and  fiber  respectively,  Gedfa  and  Ge  are  the 
gain  coefficients  of  ED  FA  and  photodetector  respectively, 
and  Ur  is  the  voltage  required  to  produce  a  phase  shift 
of  7 r  of  the  electro-optic  modulator. 


11 


Appendix  D:  Analysis  of  linear  operation  and  noise 
characteristics  of  the  optomechanical  microwave 
filter /delay  line 


In  this  section,  we  examine  the  assumption  of  weak 
mechanical  amplitude  under  strong  optical  pump  and  an¬ 
alyze  the  optomechanical  microwave  filter/delay  line  per¬ 
formance  in  terms  of  linearity  and  noise  characteristics. 
We  find  that  thermo-optic  effect  constrains  the  mechan¬ 
ical  amplitude  due  to  saturation  of  the  optomechanical 
gain.  This  effect  sets  the  linear  operation  range  and  the 
suppression  of  the  mechanical  thermal  noise. 

The  thermo-optic  effect  induced  optical  cavity  fre¬ 
quency  shift  can  be  described  by  the  following  equa¬ 
tions  [ST; 


5iuc 


/rr  sdnSi(To) 

-  wcnSi  (T0 ) - — —  AST, 


ST  = 


dT 

2 

r^c 


™Si(To)2VrPA 


(Dl) 

(D2) 


where  nc  is  cavity  photon  number,  nsi,  r,  and  ?  is  the 
refractive  index,  thermal  resistance,  and  two-photon  ab¬ 
sorption  coefficient  of  silicon  respectively,  c  is  the  speed 
of  light,  Vtpa  is  the  cavity  volume  for  two-photon  ab¬ 
sorption,  and  A  is  a  perturbation  theory  coefficient  A  = 


/si  |E(r)|2dr 


;.  Substituting  Eq. 


D2 


into  Eq. 


/ngi(To)*|E(r)|  =  dr-  __  _  _ ^ 

using  the  parameters  of  silicon  given  in  Ref.  along 
with  d  k  7.5  x  10-2,  Vtpa  ~  (VnSi (T0))3,  we  have 


Dl 


and 


Sloc  =  £n2,  £  «  -33.9  Hz. 


(D3) 


We  proceed  to  include  the  term  of  thermo-optic  effect 
(Eq.  |D3[)  to  the  equations  of  motion  (Eqs.  C3  and  C4|, 
which  are  then  modified  to  be 


iua+  =  (i(A  +  £(|ao|2  +  l«+|2)2)  -  ^)a+  ~  WoaoP-, 
-iujfi-  =  —  (iu>m  +  y)/3-  -  igoa0a\  -  y/YiPm,  An  =  V^thA 


(D4) 

(D5) 


where  we  have  denoted  a+  and  /3_  as  the  static  value  of 
the  corresponding  operators  without  input  optical  side¬ 
band  signal  and  we  have  included  the  mechanical  thermal 
noise  input.  Also,  we  specifically  consider  the  operation 
with  the  i?i  mode,  and  thus  ignore  optical  cavity  L  which 
has  much  weaker  coupling  with  R\  compared  to  optical 
cavity  R.  Eqs.|D4|and|D5|can  only  be  solved  numerically 
for  a  generic  pump  condition.  To  reveal  the  thermo-optic 
effect  on  the  mechanical  amplitude,  we  consider  a  special 
pump  condition  corresponding  to  the  original  threshold 
of  mechanical  self-oscillation,  i.e.  4gg|ao|2/(ft7,:)  =  1- 
In  this  case,  we  can  analytically  solve  for  the  down- 
converted  photon  number  and  mechanical  amplitude  at 
oj  =  wm  from  Eqs.  |D4|  and  |D5[  assuming  A  =  wm  and 
I <5+ 1  >  |<*o|, 

l«+l2  =  (^)1/5,  (D6) 

+  (D7) 

For  k  =  2tt  x  0.8  GHz,  7*  =  27r  x  3.6  MHz,  go  =  2n  x  1.39 
MHz,  nth  =  ~  1000,  we  have  |d+|2  ~  1900, 

|ao|2  ~  380,  and  |,S-|2  ~  3.7  x  105.  It  is  the  optical 
resonance  shift  induced  by  thermo-optic  effect  that  sat¬ 
urates  optomechanical  gain  and  prevents  runaway  of  the 
mechanical  amplitude  at  the  threshold. 

Now  we  can  estimate  whether  the  mechanical  ampli¬ 
tude  is  large  enough  to  induce  nonlinearity  through  exci¬ 
tation  of  higher  order  optical  sidebands.  The  nonlinear¬ 


ity  arises  due  to  pump  saturation  and  occurs  when  the 
amplitude  of  the  first  order  optical  sideband  significantly 
deviates  from  being  linearly  proportional  to  the  mechan¬ 
ical  amplitude,  i.e.  approximation  Ji(z)  ~  §  breaks 

down  [52],  where  z  =  go\ZMP-\2  +  2 /ojm  is  the  normal¬ 
ized  mechanical  amplitude.  For  |/3_|2  ~  3.7  x  105  calcu¬ 
lated  above,  the  deviation  is  only  about  1%.  In  the  ex¬ 
periment,  we  find  for  the  largest  pump  power  Ppl,r  ~  0.2 
mW,  |/3- 12  ~  2.0  x  105,  which  gives  z  =  0.18  and  a  lin¬ 
ear  deviation  of  0.4%.  As  a  result,  scattering  into  higher 
order  optical  sidebands  does  not  need  to  be  included  in 


Eq.  (D4)  and  (D5). 


Next,  we  consider  the  response  of  the  mechanical  os¬ 
cillator  to  a  small  input  optical  sideband  signal  by  per¬ 
turbative  expansion  of  Eqs.  |D4|  and  |D5|  In  this  case  the 
coherent  mechanical  amplitude  of  Eq.  |C5|  is  modified  to 
be 


P-  = 


igOLy/HeL^OiOL 


iK,-V)  +  f-£ 


(D8) 


where  5k  =  £(|aofc|2  +  \&+k\2)2  is  the  thermo-optic-effect 
induced  optical  frequency  shift.  Eqs.  |C6|  and  C8  can  be 
modified  correspondingly.  According  to  Eq.  |D8[  the  ef¬ 
fective  mechanical  loss  rate  is  7ea  =  7i  —  4g°  ^  °fe  + 


XA  49otIfc0^  (^fy^)2-  The  deviation  from  a  linear  re¬ 
sponse  can  be  caused  by  the  additional  cavity  photons 
from  the  input  signal,  and  is  characterized  by  the  ra- 
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Frequency  offset  (kHz) 

FIG.  5.  Phase  noise  suppression  ratio  between  microwave 
signal  power  of  -15  dBm  and  -30  dBm  at  largest  optical  pump 
level.  Solid  line  is  theoretical  value  calculated  using  Eq.  |D11| 
and  dots  are  experimental  data. 


tio  r  =  |a_|_fc|2/|a+fc|2  (in  our  device  the  contribution 
is  mainly  from  a+R).  At  the  theoretical  self-oscillation 
threshold,  we  find  the  1  dB  compression  point  of  the 
S— matrix  to  be  equal  to  a  microwave  power  of  -19  dBm 
(assuming  |acm|2  =  |aofl|2)-  In  the  experiment,  we  find 
for  the  largest  optical  pump  power  Ppl,r  =  0.2  mW 
(which  is  slightly  above  the  self-oscillation  threshold), 
the  1  dB  compression  point  occurs  at  a  microwave  signal 
power  of  -15  dBm.  For  reversed  operation  (cavity  R  as 
input),  the  1  dB  compression  point  is  reduced  by  a  factor 
of  70M,i/70M,i?  (assuming  |a0i|2  =  |aotf.|2)- 

We  now  analyze  the  noise  characteristics  of  the  op¬ 
tomechanical  cavity-waveguide  system.  The  dominant 
form  of  noise  is  from  thermally  excited  phonons  in  the 
system.  From  the  thermally-added  mechanical  noise  ref- 
ered  to  the  input  signal  (Ke/«;)_17inth/|70M|  [23] •  we 
define  the  noise-equivalent  optical  signal  power 

Pne  =  —  Swc— 2nB,  (D9) 

«e  70M 

where  B  is  the  bandwidth  of  the  coherent  signal  and  all 


the  quantities  are  referred  to  the  input  cavity.  Then  the 
noise-equivalent  microwave  signal  power  is 

Vne  =  A^VI  (DIO) 

7T  ip 

We  find  that  for  the  largest  pump  power  when  operating 
at  Pi  resonance,  if  the  input  port  is  cavity  L ,  the  noise 
equivalent  microwave  power  is  -30  dBm;  if  the  input  port 
is  cavity  P,  the  noise  equivalent  microwave  power  reduces 
to  -70  dBm  because  of  the  significantly  enhanced  7om  of 
cavity  P  with  Ri  mode. 

For  a  self-oscillating  mechanical  oscillator,  the  intrinsic 
oscillator  noise  can  be  suppressed  by  the  injection  of  an 
external  coherent  signal  [53] .  The  suppressed  phase  noise 
(ignoring  input  signal  noise)  can  be  modeled  by  |54) 

Sj,ct>(u)  —  U]m)  =  7eff  ,2  2‘-’W(a;  W™)>  (DU) 

1  UJ—LJm  /  P 

where  (w)  is  the  intrinsic  phase  noise  spectral  density 
without  injection  and  p  =  |/3_|/|/3_|  is  the  ratio  between 
injected  mechanical  amplitude  and  free-running  ampli¬ 
tude.  Experimentally,  we  infer  the  phase  noise  from  the 
measured  noise  power  spectral  density  using  the  defini¬ 
tion  =  Sbb(u>)  —  f  Sbb(co)dui.  At  the  largest  optical 

pump  level,  the  phase  noise  suppression  ratio  between  the 
microwave  signal  power  -15  dBm  (1  dB  S— matrix  com¬ 
pression  point)  and  -30  dBm  (noise  equivalent  power)  is 
shown  in  Fig.  [5]  The  model  of  Eq.  |Dll| explains  well  the 
measured  noise  suppression  level  in  the  offset  frequency 

—  range  between  7eff/4  and  37eff/4  (7eg  =  27r  x  17 
kHz). 

Appendix  E:  Waveguide-mediated  cavity-cavity 
coupling 

1.  Analytical  derivation 

When  7eL,ii/27r  <C  /fsr,  and  all  the  waveguide  modes 
are  outside  of  the  line  width  of  the  cavities,  two  mechan¬ 
ical  cavities  can  acquire  a  waveguide  mediated  coupling. 
The  generic  Hamiltonian  describing  this  case  is 


H  =  humLb^bL  +  hiomRVRbR  +  ^2  hukb\bk  +  ^  h(gLkb]Lbk  +  g*LkbLb{)  +  ^  HgRkbRbk  +  gRkbRb\ ),  (El) 


where  b\{bk)  is  the  creation(annihilation)  operator  of  the 
fc-th  waveguide  mode,  gLk  =  V^lJfsr  and  gRk  = 
(— )fcV/27ei?,/FSR  are  the  coupling  coefficients  of  the  left 
and  right  mechanical  cacity  modes  with  the  k— th  waveg¬ 
uide  mode,  and  the  summation  is  over  all  waveguide 
modes.  Note  the  relative  sign  between  gkk  and  gRk  comes 
from  the  symmetry  of  waveguide  modes  with  respect  to 


the  center  of  the  two  cavities. 


In  the  degenerate  cavity  case,  i.e.  =  <jJmR,  we 

can  calculate  the  coupling  V  between  the  two  cavities 
mediated  by  the  waveguide  modes  using  second  order 
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perturbation  theory  as  follows, 


V 

2tt 


E 


9Lk  9Rk 
2n  27 r 

$ ~  kfFSR 


y/'J eL^ieR 
2n 


E 


r* 


/fsr  ‘ 


—  fc7T 


y/lfeL'YeR  1 

27 r  sin-F^— 7r’ 

/FSR 


(E2) 


where  5  is  the  frequency  difference  between  the  degener¬ 
ate  cavity  modes  and  their  nearest  waveguide  mode.  For 
the  perturbation  theory  to  be  valid,  we  require 


I  gL,Rk  |/2tt  =  y/2'yeLiRfFSR/27r  «  S.  (E3) 


We  consider  two  special  cases.  First,  when  5  =  /fsr/2, 
i.e.  the  two  cavity  modes  are  in  the  center  of  two  nearest 
waveguide  modes,  then  V  =  y/'JeL'yeR-  In  the  other  limit 
when  S  <C  /fsr,  then  V  sa  ^7 eL^eR One  can  prove 

in  this  case,  V  <  v^TeL.ii /fsr  =  \gL,Rk\,  by  taking  into 
account  of  the  condition  of  Eq.  |E3[  Different  from  the 
previous  case,  here  the  contribution  to  the  coupling  can 
be  almost  exclusively  attributed  the  nearest  waveguide 
mode  [k  =  0). 


2.  Measurement  modeling 

We  model  how  to  estimate  the  phonon  waveguide  me¬ 
diated  coupling  between  two  optomechanical  cavities.  In 
the  device  for  demonstrating  waveguide-mediated  cavity 
coupling,  we  used  a  phonon  waveguide  without  air  holes. 
The  band  structure  of  the  waveguide  is  shown  in  Fig.  UK- 
Such  a  waveguide  provides  large  phonon  group  velocity 
and  thus  large  free  spectral  range  in  order  to  isolate  cav¬ 
ity  modes  from  waveguide  modes. 

In  the  sample  device  shown  in  Fig.  4  of  the  main  text, 
virtual  phonons  in  waveguide  mixes  the  mechanical  cav¬ 
ity  modes  M l  and  Mr  of  cavity  L  and  R  into  hybridized 
modes  (7+  and  C_ .  The  coupling  strength  can  be  inferred 
by  measuring  the  optomechanical  coupling  of  hybridized 
modes  C+  and  C-  with  the  two  optical  modes  Ol  and 
Or  based  on  the  following  model.  C+  and  C_  can  be  ex¬ 
pressed  as  linear  superposition  of  Ml  and  Mr,  assuming 
ignorable  energy  distribution  in  the  waveguide  (as  C+t- 
are  well  separated  from  waveguide  modes), 

C+  =  ol\Ml  +  ol^Mr,  (E4) 

C_  =p1ML+p2MR.  (E5) 

The  superposition  coefficients  cq  and  (3i  satisfy  the  fol¬ 
lowing  relation 


A,  =  ,02,  =  2|y| 

Ol  p2  &-LR  +  A_| —  ’ 


(E6) 


where  V  is  the  waveguide  mediated  coupling  between  Ml 


200  300  400  500  600  700 


FIG.  6.  a  Band  structure  of  shielded  phonon  waveguide  with¬ 
out  air  holes.  Inset  is  the  modal  profile  with  k  =  0.8n/a. 
b,c  Characterization  of  optomechanical  coupling  between  hy¬ 
bridized  mechanical  modes  C+  and  C_  with  the  optical  modes 
of  cavity  L  (b)  and  R  (c).  7  =  7i+7oM  is  effective  mechanical 
decay  rate  and  nc  is  the  cavity  photon  number. 


and  Mr,  Alr  is  the  fabrication  induced  frequency  dif- 

ference  between  Ml  and  Mr,  and  A_| _ =  yj A2LR  +  4V2 

is  the  frequency  difference  between  the  hybridized  modes 
C+  and  C-,  On  the  other  hand,  the  ratio  of  the  super¬ 
position  coefficients  is  directly  related  to  the  measurable 
optomechanical  coupling 


9+l  _  I oi,  g+R  _  ,0(2 
g-R-'fo 


(E7) 
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where,  for  example,  g+L  stands  for  the  coupling  between 
mechanical  mode  C+  and  optical  mode  Or. 

We  inferred  the  optomechanical  couplings  by  mea¬ 
suring  the  pump  dependent  effective  mechanical  damp¬ 
ing  rate  of  C+ -  (Fig.  [6Jd,c)  using  the  relation  7  = 
7 i  -  g2nc  (A-Mm)2+K2/4>  and  obtained 

g+L  =  927  ±  25  kHz,  g_L  =  710  ±  55  kHz  (E8) 
g+R  =  683  ±  90  kHz,  g_R  =  870  ±  35  kHz  (E9) 

As  a  result,  the  ratio  of  superposition  coefficients  is 

| ^-|  =  1.30  ±0.14,  |^|  =  0.79  ±0.14.  (E10) 

Pi  P2 


Along  with  A4 _ /27r  =  4.0  MHz  as  read  from  the  spec¬ 

trum,  we  find,  according  to  Eq.  |E6[ 


\V\/2ir  =  1.94  ±  0.06  MHz,  Alr/2tt  =  0.98  ±  0.48  MHz. 

(Ell) 

To  compare  with  the  analytical  formula  of  V  (Eq.  E2 1 , 


in  this  device,  we  have  /fsr  =  54  MHz  and  6  =  17 
MHz.  And  using  the  measured  V/2n  =  1.94  MHz,  we 
find  sJleLleR/2 tt  =  1.62  MHz. 


Appendix  F:  Optical  non-reciprocity  based  on 
distantly-coupled  optomechanical  cavities 

We  provide  theoretical  analysis  of  achieving  optical 
non-reciprocity  in  the  distantly-coupled  optomechanical 
cavities,  and  show  its  viability  based  on  the  typical  pa¬ 
rameters  of  our  fabricated  devices.  In  this  case,  the 
waveguide  connecting  the  two  cavities  should  support 
both  guided  mechanical  and  optical  modes,  and  the  two 
optomechanical  cavities  are  designed  to  be  identical.  As 
already  been  demonstrated  in  our  experiment  and  in 
Ref.  [2UJ  respectively,  such  a  waveguide  can  mediate  a 
tight-binding  type  coupling  for  the  mechanical  and  op¬ 
tical  cavity  modes.  The  Hamiltonian  of  the  system  can 
thus  be  written  as  follows, 


H=  ^  hujcka'kak  +  J{a'LaR  +  aLa'R)  +  E  hpmkb\bk  +  V  (b{bR  +  bLkR)  (FI) 

k—L,R  k—L,R 

+  ^2  bgk{b\  +  bk)a\ak+  ^  ih^K^apke~l“pt~l‘Pk{dk  -  a{), 
k—L.R  k=L,R 


where  J  and  V  are  the  waveguide  mediated  optical  and 
mechanical  coupling  strength,  and  the  last  two  terms  are 
the  optical  pumps  in  the  two  cavities  which  have  a  same 
frequency  and  correlated  phases. 

The  optical  non-reciprocity  arisen  from  this  system 
can  be  intuitively  understood  from  a  schematic  shown 
in  Fig.  UK  The  input  optical  signal  undergoes  a  Mach- 
Zehnder  type  of  interference  through  the  system:  one 
path  is  the  direct  photon  hopping  and  the  other  path  is 
through  radiation-pressure  interaction  induced  transition 
to  phonon  and  phonon  hopping.  The  phase  of  the  lat¬ 


ter  path  involves  the  phase  difference  of  the  two  pumps, 
which  is  tpi,  —  <pr  for  one  direction  and  (pR  —  ipi,  for  the 
reversal  direction.  Such  a  non-reciprocal  phase  resem¬ 
bles  an  effective  magnetic  flux  for  photons,  resulting  in 
the  non-reciprocal  transmission  (48j. 

We  first  consider  the  case  when  both  cavities  are  be¬ 
ing  pumped  with  blue  detuned  lasers  (uip  —  u)ck  =  uimk). 
The  equations  of  motion  of  the  system,  in  terms  of  red 
optical  sidebands  of  the  pumps,  can  be  derived  from  the 
Hamiltonian  of  Eq.  |F1|  using  rotating  wave  approxima¬ 
tion,  given  sideband  resolved  condition  iomk  Kk. 


da,L 

dt 

daR 

dt 


=  (i$L  ~  -^-)aL  ~  iJaR  -  igLaLel(l,Lb*L  -  y/K^EaL> in, 

(F2) 

=  ( idR - Y^aR  ~  ^ aL  ~  ^9RaRe^Rb*R  -  t/kIrur, in, 

(F3) 

=  -( iuimL  +  -J )bL  -  iVbR  -  igLaLelrl,La*L, 

(F4) 

=  -{iUmR  +  ~^)bR  -  iVbL  -  igROtRel(t,Ra*R, 

(F5) 

1 

where  Sk  =  ojp  —  cock  and  akel^k 

is  the  steady  state  am- 

plitude  of  the  local  optical  cavity  mode,  which  is  related 
to  the  pumping  amplitudes  as  follows 
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FIG.  7.  a  Schematic  of  the  optical  non-reciprocity  in  distantly-coupled  optomechanical  cavities  stemmed  from  the  non-reciprocal 
optical  pump  phases  and  interference  between  photonic  and  phononic  transmission  paths,  b-e  Transmission  coefficient  b(d) 
and  isolation  ratio  c(e)  for  blue(red)  detuned  optical  pumps  with  practical  device  parameters  ji  =  2n  x  3  MHz,  Ql,r  =  27t  x  0.8 
MHz,  cavity  photon  number  nCL,R  =  2000,  kl,r  =  2n  x  2  GHz,  k£l,r  =  2n  x  1  GHz,  c/>l  —  (/>r  =  ^,  J  =  2n  x  422  MHz,  and 
V  =  27 r  x  10  MHz. 


(^R(L)  —  KR(L)/^)y/^eL(R)OipL(R)e  +  iJ  ^K^^)OtpR(L)e  l<fiRW 

(iSL  -  KL/2)('iSR  -  kr/2)  +  J2 


(F6) 


We  find  the  steady  state  amplitude  is  approximately 
^/KfJ:Oipke~1,pk  /iSk  under  the  condition  <5*.  ~  0Jmk 
Kk ,  J,  which  means  each  cavity  is  effectively  only  driven 
by  its  own  optical  pump.  This  can  be  intuitively  under¬ 
stood  by  the  fact  that  even  and  odd  hybridized  optical 
cavity  modes  are  driven  equally  (as  5k  J)  and  thus  the 
amplitude  of  one  local  cavity  mode  is  cancelled  out  and 
effectively  not  being  driven.  Thus,  each  cavity-enhanced 
optomechnical  coupling  can  be  independently  controlled 
by  the  pump. 


After  solving  the  equations  of  motion,  we  calculate  the 
ratio  between  right  transmission  coefficient  TR  and  left 
transmission  coefficient  TR  of  the  optical  signal  and  find 


Tr 

Tl 


J  - 
J  - 


_ _ 

(i(u>-uirnL)+:^-)(i(u>-ujrrl,R)+^-)+V2 


it<t>L  —  4>r) 


VgLgR&LttR 


g  —  i(4>L  —  <t>R) 


(F7) 

Interestingly,  this  ratio  is  not  explicitly  dependent  on  5k 
and  Kk  as  an  intrinsic  property  of  the  device.  At  the 
poles  to  =  ( uimL  +  umR  ±  y/ ( wmL  -  UrnR.)2  +  Y2)/2,  he. 
frequency  of  the  hybridized  mechanical  modes,  and  as¬ 


suming  t  7i,  we  have 

T  4-  _ VgLgnaLOLR _ 0i(<j>L-^R.) 

Tr  = 

Tl  j  =f  i _  VgMMMm  c-i(<b r  -^.r 

'Hy/V2+{UmL-UmR)2  /4 

Thus  perfect  non-reciprocity,  i.e.  one  direction  has  van¬ 
ishing  transmission  while  the  other  direction  has  maximal 
transmission,  can  be  achieved  by  satisfying  the  following 
condition 


4>l  —  <t>R  =  ± 


7T 

2’ 


J  = 


_ VgLgRQLaR _ 

li\/V2  +  ( u}mL  -  ujm_ r)2/4 


(F9) 

Under  this  condition,  the  transmission  coefficient  for  the 
through  direction  is  (for  simplicity  assuming  c dmL  — 
M mR ) 


^7^0  —  yj  KeL^eR 


2  9l_9r^l^R. 

7  i 


(  Kl  _  9 L  9 R®1- L  CX- R 
V  2  7  i 


KR  _  9 L  g R&LOlr 
2  7?- 


)' 

(F10) 

Similar  results  can  be  derived  for  the  case  of  red  de¬ 
tuned  pumps,  and  we  find  the  perfect  non-reciprocity 
condition  (Eq.  F9 1  remains  the  same;  while  the  trans¬ 
mission  coefficient  for  the  through  direction  at  poles  is 
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given  by 


T^o  =  y/^eL^eR 


2  9l9r^l^R. 
7  * 


(*l 
V  2 


9 L  9 R& L& Ft  \(  KR 
'Vo  /  v  2 


9l9r<^l(^r 

7  i 


)' 

.  (Fn) 

We  note,  in  general,  an  amplified  transmission  in  blue 
detuned  case  and  an  attenuated  transmission  in  red  de¬ 
tuned  case  for  the  through  direction  at  poles.  One 
can  prove  from  Eq.  |F11|  that  in  the  red  detuned  case, 
T^o  <  y  KeL^eR./ {klKr)  <  1  and  equality  is  achieved 
when  Kk /2  =  Srss^lor  .  Comparing  to  Eq.  |F9l  the  max¬ 
imal  transmission  efficiency  is  achieved  when  loss  rate 
Kk/ 2  and  coupling  rate  J  is  matched  at  the  two  cavities 
(note  we  used  umi  =  u!mf j  for  Eqs.  F10  and  Fll  I. 


Based  on  our  fabricated  devices,  realizing  the  condi¬ 
tions  for  perfect  non-reciprocity  (Eq.|F9[)  is  quite  promis¬ 


ing.  For  a  typical  device  with  7*  =  2-7T  x  3  MHz,  g l  = 
gn  =  2n  x  0.8  MHz,  maximal  available  cavity  photon 
number  nCL  =  nCR  =  2000,  and  assuming  uJmL  =  uimn, 
Eq.  |F9|  determines  Jmax  =  422  MHz.  Thus,  as  long  as 
J  <  Jmax  in  this  device,  perfect  non-reciprocity  condi¬ 
tion  can  always  be  achieved  by  tuning  the  pump  power 
and  phase.  We  have  numerically  simulated  waveguide 
(without  acoustic  shielding)  mediated  optical  coupling 
between  two  optomechanical  cavities  to  be  500  MHz  and 
less.  Thus  it  is  indeed  viable  to  demonstrate  optical  non¬ 
reciprocity  in  our  devices.  Using  these  parameters  along 
with  5k  =  L*Jmk,  Kk  =  27r  x  2  GHz,  Kek  =  2-7T  x  1  GHz  and 
<f>L  -  <i>R  =  § ,  j  =  2tt  x  422  MHz,  V  =  2tt  x  10  MHz, 
we  plotted  the  transmission  coefficients  and  isolation  ra¬ 
tio  as  calculated  from  Eqs.  |F2|fF5|  for  the  blue-detuned 
pump  In  Fig.  [7}o,c  and  red-detuned  pump  In  Fig.  0he, 
respectively. 


